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' A qualitative analysis of the chiral phase transition in QCD with two massless quarks and non- 

zero baryon density is performed. It is assumed that at zero baryonic density, p = 0, the temperature 
phase transition is of the second order and quark condesate 77 = (0 | itit | 0) = (0 | dd | 0) may be 
taken as order parameter of phase transition. The baryon masses strongly violate chiral symmetry, 

■ nriB ~ (0 I I 0)^^^. By supposing, that such specific dependence of baryon masses on quark 
^ ' condensate takes place up to phase transition point, it is shown, that at finite baryon density p the 

QO , phase transition becomes of the first order at the temperature T — T'ph(p) for p > 0. At temperatures 

CO ' T'cont(p) > T > Tph(p) there is a mixed phase consisting of the quark phase (stable) and the hadron 

' phase (unstable). At the temperature T — Tcont{p) the system experiences a continuous transition 

to the pure chirally symmetric phase. 

O . 
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^ ■ 1 Introduction 

Oh. It is well known, that the chiral symmetry is valid in perturbative quantum chromodynamics (QCD) 



X 



with massless quarks. It is expected also, that the chiral symmetry takes place in fuU-perturbative 
and nonperturbative QCD at high temperatures, (T > 200 MeV), if heavy quarks {c,b,t) are ignored. 
The chiral symmetry is strongly violated, however, in hadronic matter, i.e. in QCD at T = and low 
density. What is the order of phase transition between two phases of QCD with broken and restored 
' chiral symmetry at variation of temperature and density is not completely clear now. There are different 
opinions about this subject (for a detailed review see Ref. and references therein). 

In this talk I discuss the phase transitions in QCD with two massless quarks, u and d. Many lattice 
calculations [21 12 El El indicate, that at zero chemical potential the phase transition is of the second 
order. It will be shown below, that the account of baryon density drastically changes the situation 
and the transition becomes of the first order, and, at high density, the matter is always in the chirally 
symmetric phase. 

The masses of light u, d, s quarks which enter the QCD Lagrangian, especially the masses of u and 
d quarks, from which the usual (nonstrange) hadrons are built, are very small, mu,md < 10 MeV as 
compared with characteristic mass scale M ^ 1 GeV. Since in QCD the quark interaction proceeds 
through the exchange of vector gluonic field, then, if light quark masses are neglected, QCD Lagrangian 
(its light quark part) is chirally symmetric, i.e. not only vector, but also axial currents are conserved and 
the left and right chirality quark fields are conserving separately. This chiral symmetry is not realized 
in hadronic matter, in the spectrum of hadrons and their low energy interactions. Indeed, in chirally 
symmetrical theory the fermion states must be either massless or degenerate in parity. It is evident, 
that the baryons (particularly, the nucleon) do not possess such properties. This means, that the chiral 
symmetry of QCD Lagrangian is not realized on the spectrum of physical states and is spontaneously 
broken. According to Goldstone theorem spontaneous breaking of symmetry leads to appearance of 
massless particles in the spectrum of physical states - the Goldstone bosons. 

Let us first consider a case of the zero baryonic density and suppose that the phase transition from 
chirality violating phase to the chirality conserving one is of the second order. The second order phase 
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transition is, generally, characterized by the order parameter rj. The order parameter is a thermal average 
of some operator which may be chosen in various ways. The physical results are independent on the choice 
of the order parameter. In QCD the quark condensate, rj = |(0|mu|0)| = |(0|(i(i|0)| > 0, may be taken 
as such parameter. In the phase of broken chiral symmetry (hadronic phase) the quark condensate is 
non-zero and has the normal hadronic scale, (0 | | 0} = {25QMeV)^ , in the phase of restored chiral 
symmetry it is vanishing. 

The quark condensate has the desired properties: as it was demonstrated in the chiral effective 
theory 13 |H1, decreases with the temperature increasing and an extrapolation of the curve ri(T) to 
the higher temperatures indicates, that 77 vanishes at T = Ti°^ 180 MeV. Here the superscript "0" 
indicates that the critical temperature is taken at zero baryon density. The same conclusion follows from 
the lattice calculations E], where it was also found that the chiral condensate r] decreases with 
increase of the chemical potential |1U[ [TT] . 

Apply the general theory of the second order phase transitions (12) and consider the thermodynamical 
potential ^(jj) at the temperature T near Tc^\ Since rj is small in this domain, ^(rj) may be expanded 
in 77: 

*(^) -*o + ^^?7' + Jsr;4, B>0. (1) 

For a moment we neglect possible derivative terms in the potential. 

The terms, proportional to rj and 77^ vanish for general reasons (T2|. In QCD with massless quarks the 
absence of rj and rj'^ terms can be proved for any perturbative Feynman diagrams. At small t = T — Tc^'^ 
the function A{t) is linear in t: A{t) = at, a > 0. If i < the thermodynamical potential $(77) is minimal 
at 7; 7^ 0, while aX t > the chiral condensate vanishes 77 = 0. At small t the t-dependence of the 
coefficient B{t) is inessential and may be neglected. The minimum, fj, of the thermodynamical potential 
can be found from the condition, d^/drj = 0: 

r v/^W^, <<0; , . 

0, t>0. 

It corresponds to the second order phase transition since the potential is quartic in rj and - if the derivative 
terms are included in the expansion - the correlation length becomes infinite at T = Ti*^-* . 



2 Nucleon mass and quark condensate 

I show now, that the existence of large baryon masses and the appearance of violating chiral symmetry 
quark condensate are deeply interconnected and even more, that baryon masses arise just due to quark 
condensate. I will use the QCD sum rule method invented by Shifman, Vainstein and Zakharov [13], in 
its applications to baryons [14-17]. (For a review and collection of relevant original papers see [18]). The 
idea of the method is that at virtualities of order ~ 1 GeV^ the operator product expansion (OPE) 
may be used in consideration of hadronic vacuum correlators. In OPE the nonoperturbative effects reduce 
to appearance of vacuum condensates and condensates of the lowest dimension play the most important 
role. The perturbative terms are moderate and do not change the results in essential way, especially in 
the cases of chiral symmetry violation, where they can appear as corrections only. 

For definiteness consider the proton mass calculation [14-16]. Introduce the polarization operator 

n(p) = 7 J d'xe'p^o\Tax),mm m 

where ^(x) is the quark current with proton quantum numbers and is chosen to be space-like, 

< 0, Ip^I 1 GeV^. The current ^ is the colourless product of three quark fields, ^ — e"''^ q'^q^q'^, q = 
u, d, the form of the current will be specialized below. The general structure of n(p) is 

n(p)=p/i(p) + /2(p) (4) 

The first structure, proportional to p is conserving chirality, while the second is chirality violating. 
For each of the functions fi{p'^), i — 1,2 the OPE can be written as: 
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Figure 1. The bare loop diagram, contributing 
to chirality conserving function /i(p^): solid lines 
correspond to quark propagators, crosses mean 
the interaction with external currents. 



Figure 2. The diagram, corresponding to chi- 
rality violating dimension 3 operator (quark con- 
densate). The dots, surrounded by circle mean 
quarks in the condensate phase. All other nota- 
tion is the same as on Fig.l. 



fiip') = Yl C«(P')(0|0«|0) (5) 

n 

where (0|oi'^|0) are vacuum expectation values (v.e.v) of various operators (vacuum condensates), C'n^ 
are coefficient functions calculated in QCD. For the first, conserving chirality structure function fi{p'^) 
OPE starts from dimension zero {d = 0) unit operator. Its contribution is described by the diagram of 
Fig.l and 

pMp^) = CoPp'^ln[Al/{-p'^)] + polynomial, (6) 

where Co is a constant, A^^ is the ultraviolet cutoff. The OPE for chirality violating structure f2{p^) 
starts from d = 3 operator, and its contribution is represented by the diagram of Fig. 2: 

/2(p2) = Cip'^^o\Oqq\o)in-^-^ + polynomial (7) 

Let us for a moment restrict ourselves to this first order terms of OPE and neglect higher order terms 
(as well the perturbative corrections). 

On the other hand, the polarizaion operator (3) may be expressed via the characteristics of physical 
states using the dispersion relations 

1 f Imfi(s') 

fi{s) = — / — — ^ ds' + polynomial, s = —p^ (8) 

IT J s + s 

The proton contribution to /mn(p) is equal to 

Imn{p) = 7r(0|eb)(p|?|0)5(p2 - m^) = TrX%{p + m)6{p^ - m^), (9) 

where 

m\p) = xnv{p), (10) 

Aat is a constant, v{p) is the proton spinor and m is the proton mass. Still restricting ourselves to this 
rough approximation, we may take equal the calculated in QCD expression for Il{p) (Eq.'s(6),(7)) to its 
phenomenological representation Eq.(9). The best way to get rid of unknown polynomial, is to apply to 
both sides of the equality the Borel (Laplace) transformation, defined as 



■*oo ,s—*oo, 
s/n = M2 = Co 



)n oo 
n 



BM^m= ^ -rf^l /(«) = r / dslmfis)e-^f^' (11) 





if /(s) is given by dispersion relation (8). Notice, that 
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Owing to the factor l/(n — 1)! in (12) tlic Borcl transformation suppresses the contributions of high 
order terms in OPE. 

Specify now tlie quark current ^(x). It is clear from (9) that proton contribution will dominate in 
some region of the Borel parameter ^ rn? only in the case when both calculated in QCD functions /i 
and /2 are of the same order. This recjuirement, together with the requirements of absence of derivatives 
and of renormcovariance fixes the form of current in unique way (for more details see [14,16]): 

^{x) = e«^^(w«C7^u^)7^75rf" (13) 

where C is tlie charge conjugation matrix. With the current t/(.x) (13) the calculations of the diagrams 
Fig. 2 can be easily performed, the constants Co and Ci are determined and after Borel transformation two 
equations (sum rules) arise (on the phenomenoloical sides of the sum rules only proton state is accounted) 

M6 = A^^e-™'/^' (14) 
- 2{2Trf{0\qq\0)M^ = m~X% e"™'/*^' (15) 
= 327r* A^ 

It can be shown that this rough approximation is valid at M w ro. Using this value of M and dividing 
(15) on (14) we get a simple formula for proton mass [14]: 

m= [-2(27r)2(0|gg|0)]V3 (16) 

This formula demonstrates the fundamental fact, that the appearance of the proton mass is caused by 
spontaneous violation of chiral invariance: the presence of quark condensate. (Numerically, (16) gives 
the experimental value of proton mass with an accuracy better than 10%). 

A more refined treatment of the problem of the proton mass calculation was performed: high order 
terms of OPE were accounted, as well as excited states in the phenomenological sides of the sum rules 
and the stability of the Borel mass dependence was checked. In the same way, the hyperons, isobar and 
some resonances masses were calculated, all in a good agreement with experiment [14,17]. 1 will not dwell 
on these results. The main conclusion is: the origin of baryon masses is in spontaneous violation of chiral 
invariance - the existence of quark condensate in QCD. 

3 The chiral phase transition in the presense of finite 
baryon density 

I would like to consider here the influence of baryon density on the chiral phase transition in hadronic 
matter. Kogan, Kovner and Tekin [19] have suggested the idea, that baryons may initiate the restoration 
of chiral symmetry, if their density is high - when roughly half of the volume is occupied by baryons. 
The physical argument in favour of this idea comes from the hypothesis (supported by calculation in 
chiral soliton model of nucleon [20]), that inside the baryon the chiral condensate has the sign opposite 
to that in vacuum. This hypothesis is not proved. Even more, it is doubtful, that the concept of quark 
condensate inside the nucleon can be formulated in a correct way in quantum theory. But the idea on the 
strong influence of baryon density on the chiral phase transition looks very attractive. For this reason 
no assumption on the driving mechanism of chiral phase transition at zero baryon density will be done 
here. The problem, under consideration is: how the phase transition changes in the presence of baryons. 
The content of my talk closely follows ref.21. 

For completeness it must be mentioned, that the variation of quark codensate with temperature is 
not the only source of baryon effective mass shift. At low T the effective baryon mass shift arises also 
due to interaction with pious in thermal bath [22,23]. However, this mass shift, which may be called 
external (unlike the internal, arising from variation of quark condensate) is related to effective mass, i.e. 
to propagation of baryon in the matter and has nothing to do with the properties of the matter as a whole 
and the phase transition. (A similar phenomenon takes place in case of vector mesons, where because of 
interaction with pions in the thermal bath the mixing with axial mesons arises [24].) 



4 



Consider the case of the finite, but small baryon density p (by p we mean here the sum of baryon 
and anti-baryon densities). For a moment, consider only one type of baryons, i.e. the nucleon. The 
temperature of the phase transition, Tph, is, in general, dependent on the baryon density, Tph = Tph(p), 
with Tph(yO = 0) = ri"^ At T < Tph(/3) the term, proportional to Ep, where E = ^ + m? is the 
baryon energy, must be added to the thermodynamical potential As was shown above the nucleon 
mass m (as well as the masses of other baryons) arises due to the spontaneous violation of the chiral 
symmetry and is approximately proportional to the cubic root of the quark condensate: m = crj^^^, with 
c = (Stt^)^/^ for a nucleon. At small temperatures T the baryon contribution to $ is strongly suppressed 
by the Boltzmann factor e~^/'^ and is negligible. Below we assume that the proportionality m ~ r]^/^ 
is valid in a broad temperature interval. Arguments in favor of such an assumption are based on the 
expectation that the baryon masses vanish at T = ?ph(/o) and on the dimensional grounds. Near the 
phase transition point E — -^/p^ + m? « p + if/^ /(2p). At 77 — > all baryons are accumulating near 
zero mass and a summation over all baryons gives us - instead of eq. - the following: 

$(77, p) - $0 + \at + is ^4 ^ Cjf'^p , (17) 

where C — '^icf/{2pi). The term p^^Pi is absorbed into $0 since it is independent on the chiral 
condensate ry. The typical momenta are of the order of the temperature, pi ~ T . Thus, Eq. (|17|l is valid 
in the region r] <^T^ . In the leading approximation the coefficient C can be considered as independent 
on the temperature at T ~ TjP^ . 

Due to the last term in Eq. (|17|l the thermodynamical potential always have a local minimum at 77 = 
since the condensate 77 is always non-negative. At small t < there also exists a local minimum at 77 > 0, 
which is a solution of the equation: 

^^(^at + Brf)^+lcpr"'^Q- (18) 




Figure 3: (a) Graphical representation of Eq. p8|l : "I" is the first term and "H" the second term (with 
the opposite sign) in the r.h.s. of the equation, (b) The thermodynamic potential H17|l vs. the chiral 
condensate at a fixed baryon density p > 0. At low enough temperatures, T = Ti, the system resides 
in the chirally broken (hadron) phase. The first order phase transition to the quark phase takes place 
at Tph = T2 > Ti. At somewhat higher temperatures, T3 > Tph the system is in a mixed state. The 
temperature T4 = Tcont corresponds to a continuous transition to the pure quark phase, in which the 
thermodynamic potential has the form T5. 

At small enough baryon density p, Eq. (|18|l [visualized in Figure Ol^a)] has, in general, two roots, 
771 < TjQ and 772 > 770, where 770 = {—at/iB)^/^ is the mininmm of the first term in the right-hand side 
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of Eq. H18|) . The calculation of the second derivative d^^/drj^ shows that the second root 772 (if exists) 
corresponds to minimum of and, therefore, is a local minimum of $. The point rj — rji corresponds 
to a local maximum of the thermodynamical potential since at this point the second derivative is always 
non-positive. 

The thermodynamical potential $(77, p) at (fixed) non-zero baryon density p has the form plotted in 
Figure |3|[b) . At low enough temperatures (curve Ti) the potential has a global minimum at > and 
system resides in the chirally broken (hadron) phase. As temperature increases the minima at = and 
at 77 = 7^2 > becomes of equal height (curve T2 = Tph). At this point the first order phase transition to 
the quark phase takes place. At somewhat higher temperatures, T — > Tph, the 77 > minimum of 
the potential still exist but $(77 = 0) < $(^2)- This is a mixed phase, in which the bubbles of the hadron 
phase may still exist. However, as temperature increases further, the second minimum disappears (curve 
T4 = Tcont). This temperature corresponds to a continuous transition to the pure quark phase, in which 
the thermodynamic potential has the form T^. 

Let us calculate the temperature of the phase transition, rpii(p), at non-zero baryon density p. The 
transition corresponds to the curve T2 in Figure |2Ib), which is defined by the equation $(7^2, p) ~ ^{rj = 
0, p), where 772 is the second root of Eq. as discussed above. The solution is 

/ \ 3/5 / X 2/5 

-.'<^)-i"'-^(^) (!) . (-) 

1/2 

and the second minimum of the thermodynamic potential is at 772 = [4a (T'"-* — Tph(p))/(5 B)] 

At a temperature slightly higher than Tph(p) the potential is minimal at 77 = 0, but it has also an 
unstable minimum at some 77 > 0. The existence of metastable state is also a common feature of the first 
order phase transition {e.g., the overheated hquid in case of hquid-gas system). With a further increase 
of the density p (at a given temperature) the intersection of the two curves in Figure 3(a) disappears and 
the two curves only touch one another at one point rj = 7)4. At this temperature a continuous transition 
(crossover) takes place. The corresponding potential has the characteristic form denoted as T4 in Figure 
3(a). The temperature T4 = Tcont is defined by the condition that the first (|18(l and the second derivatives 
of Eq. (fTTjl vanish: 

/ \ 3/5 / X 2/5 

and the value of the chiral condensate, where the second local minimum of the potential disappears is 

given by 774 = [2a(rcont(p) ~ 2c )/(5-B)] . At temperatures T > Tcont (p) the potential has only one 
minimum and the matter is in the state with the restored chiral symmetry. Thus, in QCD with massless 
quarks the type of phase transition with the restoration of the chiral symmetry strongly depends on the 
value of baryonic density p. At a fixed temperature, T < Tc'^\ the phase transition happens at a certain 
critical density, pph. According to Eq. H19|l the critical density has a kind of a "universal" dependence 
on the temperature, pph(T) oc [Tc^'' — T]^^^, the power of which does not depend on the parameters of 
the thermodynamic potential, a and B. 

The expected phase diagram is shown qualitatively in Figure 4(a). This diagram does not contain 
an end-point which was found in lattice simulations of the QCD with a finite chemical potential [25,26]. 
One may expect that this happens because in our approach a possible influence of the confinement on the 
order of the chiral restoration transition was ignored. Intuitively, it seems that at low baryon densities 
such influence is absent indeed: the deconfinement phenomenon refers to the large quark-anti-quark 
separations while the restoration of the of the chiral symmetry appears due to fluctuations of the giuonic 
fields in the vicinity of the quark. However, the confinement phenomenon dictates the value of the baryon 
size which can not be ignored at high baryon densities, when the baryons are overlapping. If the melting 
of the baryons happens in the hadron phase depicted in Figure 4(a), then at high enough density the 
nature of the transition could be changed. This may give rise in appearance of the end-point observed in 
Ref. [25,26]. 

The domain where the inequality \at\ 3> Cprj^^^, p 7^ is fulfilled, has specific features. In this 
domain the phase transition looks like a smeared second order phase transition: the specific heat has 
(approximately) a discontinuity at the phase transition point, ACp = a^Tc/2B. This statment follows 
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Figure 4: The qualitative phase diagram at finite baryon density and temperature based on the analysis 
(a) without and (b) with indication of the approximate 2-nd order transition domain. 



from general theory [12]. At | at |3> Cpif'/^ the last term in (17) may be neglected and we find for the 
entropy 

(21) 



IdA ^ 



In the phase above phase transition, rj — 0, S — Sq. Below phase transition 

S = So + ^{T-T,) 

The specific heat Cp = T[dS/dT)p below the phase transition in the hmit T —> Tc is equal 

^2 



Cp — Cp„ + 2j^^'^ 



(22) 



(23) 



The correlation length increases as 

(J _ T^">) " at T - Ti°^ -> 0. The latter arises if we include 
the derivative terms in the effective thermodynamical potential. The phase diagram with this domain 
indicated may look as it is shown in Figure 4(b). Note that the applicability of our considerations is 
limited to the region |T - ri^V^'i"^ < 1 and low baryon densities. 

In the real QCD the massive heavy quarks (the quarks c, 5, t) do not influence on this conclusion, since 



their concentration in the vicinity of T « 



(0) 



200 MeV is small. However, the strange quarks, the 



mass of which mj « 150 MeV is just of order of expected t!P\ may change the situation. This problem 
deserves further investigation. 

This work was supported in part by INTAS grant 2000-587, RFBR grants 03-02-16209 and by funds 
from EC to the project "Study of Strongly Interacting Matter" under contract 2004 No. R113-CT-2004- 
506078. 
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